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We reexamine the functional renormalization-group theory of wetting transitions. 
As a starting point of the analysis we apply an exact equation describing renormal- 
ization group flow of the generating functional for irreducible vertex functions. We 
show how the standard nonlinear renormalization group theory of wetting transi- 
tions can be recovered by a very simple truncation of the exact flow equation. The 
derivation makes all the involved approximations transparent and demonstrates the 
applicability of the approach in any spatial dimension d > 2. Exploiting the non- 
uniqueness of the renormalization-group cutoff scheme, we find however, that the 
capillary parameter a; is a scheme-dependent quantity below d = 3. For d = 3 the 
parameter uj is perfectly robust against scheme variation. 
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INTRODUCTION 

Wetting transitions [l, 2] have posed a considerable challenge to condensed matter theory 
over the last decades. Of particular focus is the controversial critical wetting transition [3] 
in systems with short-ranged intermolecular interactions. This is of remarkable interest to 
the renormalization-group (RG) theory. The physical dimensionality d = 3 turns out to be 
the upper critical dimension yielding subtle, non-universal behavior very distinct to bulk 
criticality. The transitions in question are specific in a number of respects: (1) critical 
exponents are non-universal; (2) fluctuation effects are important in d £ [2, 3] even though 
anomalous scaling of the corresponding propagator is absent (the critical exponent r] = 0); 
(3) even though fluctuation effects are important, they give rise to no phase boundary shifts 
in many physical cases in d close to or equal three; (4) in the so-called strong-fluctuation 
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regime, the interfacial correlation length diverges exponentially at the transition in d = 3 (at 
the upper critical dimension!), while for d G [2,3[ and above d = 3 it displays a power-law 
singularity. 

Renormalization group theory has proven an efficient tool to analyze the remarkable 
fluctuation effects near the wetting transitions. The early studies 4], l5| were performed in 
different formulations and included only terms linear in the effective interfacial potential. 
Although such an approach ultimately yields correct critical behavior in d = 3, it fails to 
controllably describe the strongly repelling interaction between the interface dividing the 
coexisting phases and the solid wall at short separations. For d < 3 nontrivial fixed points 
describing the disordered phase and the critical manifold in the functional space are missed 
by a linearized RG framework. So are fixed points describing multicritical behavior. 

A significant improvement was offered in Refs. 0, Q], where a non-linear RG procedure 
was carried out relying on Wilson's approximate recursion relations. The study yields a 
general, unified description of a wide class of models featuring wetting-type transitions. 
In particular, it shows that the transitions in the strong-fluctuation regime, occurring for 
quite distinct interfacial models, are described by just one fixed point. The study was 
subsequently extended to account for multicritical phenomena [8|. This theory relies on 
approximate recursion relations, derived following Wilson's scheme making a number of 
approximations. Below we rederive the continuous version of these relations using the one- 
particle irreducible variant of functional RG. The derivation is very simple and relies on clear 
and modest assumptions. We subsequently exploit the ambiguity in defining the RG scheme. 
We analyze the linear RG limit and argue that the capillary parameter u, influencing critical 
properties and determining the boundaries between different scaling regimes is a scheme- 
dependent quantity below d — 3. Subsequently we show that no such dependence occurs at 
the physical dimensionality d = 3. 

THE EXACT RG FLOW EQUATION 

In the context of critical wetting phenomena, we consider a capillary-wave type effective 
Hamiltonian 

(1) 



H[l] = H\l}/k B T = j d d ~ x x °- (VZ~) 2 + Vil) 
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where Z(x) describes the position of a fluctuating interface dividing two coexisting phases 
above a point x on a planar, uniform wall. The interfacial stiffness a controls the energetic 
cost of deforming the interface with respect to a planar configuration. The effective interfa- 
cial potential V(l) is infinite or very large for / < to prevent the interface from penetrating 
the wall. It displays a minimum at / = Z corresponding to the mean-field (MF) value of the 
interfacial separation, and approaches a constant value as I — > oo. By varying temperature 
or another control parameter the system may be tuned to a wetting transition, where the 
average value of I continuously diverges. In the MF and weak-fluctuation scaling regimes 
(see p, Hi) this coincides with the vanishing of the attractive tail in Vil). The capillary-wave 
models are capable of describing only large-scale interfacial physics. The variables should be 
understood as coarse-grained over a scale of the order of the intrinsic width of the interface, 
which in turn is of the order of the bulk correlation lengths. The Hamiltonian Eq. ([1]) is 
therefore supplemented by an upper momentum cutoff A . 

The theoretical framework we employ is usually referred to as one-particle irreducible 
functional RG, or effective action RG. For the remaining part of the section we outline a 
derivation of an exact flow equation, following Refs. . 

We first consider the functional 

w A [h] = log J V i e -m-AH A ]i]+jd^ X h^m (2) 

When A%a[Z] vanishes, Wa [h] is just the generating functional for the connected correlation 
functions. The additional term AHa[1] is assumed quadratic in I 

A^A0 = \f d d - l x J d^yR^yYl^liy) = l - J -|^i? A ( g )f q L q . (3) 

The cutoff function R\(q) is required to vanish when the cutoff scale A is sent to zero, and 
become infinite (or very large) for A — > Aq at any fixed wavevector q. Due to the presence 
of the term AH a, the Fourier modes of I corresponding to sufficiently small q (q < A) aquire 
an effective g-dependent mass providing an infra-red cutoff. The short modes (q > A) are 
not affected by the term AT-La- 

The scale-dependent effective action is defined via the modified Legendre transform: 

r A [/(x)] = -W A [h{x)] + J d d ~ l xh( x )l( x ) - AH A [l(x.)} , (4) 



where 



<M = <r<*» = ^. w 
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The last term in Eq. PJ provides that T\ [l] = (see e.g. Ref. |10(] ) . The quantity T\[l] 
interpolates between the Hamiltonian (for A = A ) and the usual Legendre transform of the 
logarithm of the partition function, i.e. the Gibbs free energy for A — > 0, when the infrared 
cutoff is removed. Executing the derivative of Eq. (SJ) with respect to A, after simple algebra 



we obtain the following flow equation [9| 

arm- 1 f "til 3 a-Ra(9) , 

where 

Solution of Eq. ([6]) with the Hamiltonian T-L[l] as the initial condition is equivalent to eval- 
uating the corresponding partition function. 

TRUNCATION 

The functional differential equation Eq. can hardly ever be solved exactly. However, 
in many contexts it provides an interesting starting point for non-standard approximation 
schemes that may go beyond perturbation theory. Below we show how the basic equation of 
the non-linear theory of wetting is simply recovered from Eq. (J6]) relying on two assumptions: 
(1) locality of the effective action, (2) neglecting the renormalization of the stiffness coeffi- 
cient & and higher order gradient terms. The discussion of these assumptions is postponed 
to the last section. 

Relying on the assumed locality, we expand the effective action in gradients of /(x). 

1, 



r A [z] = J d d ~ l x 



U A (l) + ^S A (/)(V/) 2 + .. 



(8) 



and subsequently neglect the renormalization of the gradient terms, putting £a(0 = a f° r 
all A. Fully analogous approach was developed for bulk criticality described by the O(N) 



models 1CH12|. In that context the order parameter field shows an anomalous dimension 
at criticality (the exponent 77 7^ 0) for d < 4 and therefore the renormalization of the 
gradient term should not be neglected. Accurate calculation of the critical exponents with 
this method in d = 3 requires accounting for terms of the order <9 4 . The present situation 
is completely different, as r\ = for any d at the wetting transition [3]. What is neglected 
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is therefore only a finite renormalization of the stiffness coefficient. The approach should 
therefore yield very accurate results. 

By plugging the derivative expansion Eq. (JEJ) into Eq. (jBJ) we arrive at a flow equation for 
the effective potential U\(l): 

OaU a {1) - 2 J — 2 - um - RA(q) , (9) 

where 

Kd) = *™ . do) 

We now introduce the following variables: 

» = -iogp i/ = £. * = ^a^, «a(«) = tr A (i)A- w , r A ( y ) = ^M, (ii) 

A A z oq £ 

in terms of which Eq. (Q takes the scale-invariant form: 

dsM*) = J dyy^ 1 + + rA(y ) +(d-lM*)+^*^A(*) ' (12) 

d+1 

where we performed integration over the angular variables, and Sd = ^fxpk is t he surface area 
of a (^-dimensional sphere of radius one. We stress, that Eq. ( fl2l) was derived assuming not 
more than that a remains constant during the flow and the local structure of the effective 
action is pre served by the flow. The choice of the cutoff function is to a large extent 



arbitrary 



10] . As long as no approximations are made, the results obtained in the limit 



A — > do not depend on the particular form of R A (q). Truncation of the exact flow 
equation introduces a weak cutoff dependence of the computed quantities, which can serve 
the purpose of estimating the errorbars. 

We now specify the limit, in which the standard theory is recovered. We consider the 
sharp momentum cutoff, where r\(y) — >■ for y > 1 and r\(y) — > oo for y < 1. This is 
realized by taking 

R A (q) = lim ctAV(A 2 - q 2 ) ■ (13) 

7— +00 



With this choice, the remaining integral in Eq. ([12]) can be done. This yields 

S 3 c/ 
d s u K (z) = 2(2n)d-i + 9 l u ^( z )) + (d~ ^) u k{z) H —zd z u A (z) + const , (14) 

where the constant term shifts the flowing potential u A {z) without influencing its shape and 
can therefore be skipped. After making the following substitution of variables: 



MA= ;^p' ^V^p 2 '' (15) 
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we obtain 



K'^bT , M 1 ^ / m , ^ / a , 3-d , 



M*0 = (47r)( ;_ 1)/2r( ^ 1) Mi + ^xf >(*')) + (d - iM*0 + — s^CO > (16) 

which coincides with the continuous limit of the Wilson's approximate recursion relations 

_,-,„. .-a 

we recover the linear flow equation analyzed in Ref. [5[. 

Although the present analysis yields the same result as Ref. |t| it avoids approximations 
made while performing the partial trace over the short modes in the Wilson scheme. The 
"arbitrary" length scale, which in [3| is determined using phenomenological arguments and 
so that the linear RG [5] is recovered in the appropriate limit is absent here. We evaded 
the need for any approximations regarding the eigenfunctions into which the short modes 
are expanded. In the present framework, the basic nonlinear RG equation of wetting phe- 
nomena turns out to emerge in the sharp-cutoff limit from an exact flow equation merely by 
demanding that the stiffness coefficient is not renormalized. Also note, that the linear RG 
equation [5| is recovered by neglecting terms only of the order (d^,v(z')) 2 . In fact, the order 
of these terms does not depend on the choice of the cutoff. This yields the conclusion, that 
the reliability of the linear RG approach does not depend on the magnitude of V(l) and 
djV{l) in the vicinity of the wall, and the neglected terms are only quadratic in curvature. 



CUTOFF SCHEME DEPENDENCIES 



We now pass to another point of this paper. It has to be emphasized, that despite 
the modesty of the approximations made, the equation defining the non-linear theory of 
wetting cannot be regarded as unique. Indeed, the Lipowsky-Fisher framework is recovered 
only in the sharp cutoff limit given by Eq. ( 1X3]) . Exploiting the freedom of the choice of 
the cutoff function R\(q), a whole family of flow equations can be generated. Variation of 
the scale parameter A corresponds to integrating out a different set of modes depending 
on the cutoff R\(q). The fact that different regulators yield distinct flow equations is not 
particularly surprising since for two distinct choices of Ra{q) the corresponding effective 
potentials considered at a scale A describe two different physical systems. The present 
formulation of functional RG assures that the RG trajectories are universal in the so called 
freezing regime, once the cutoff scale is lowered below the scale set by renormalized masses 
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present in the problem and the tree contributions to the beta-functions dominate. For a 
detailed discussion as well as connection to other exact RG frameworks (Wilson-Polchinski 



and Wegner- Houghton) see Ref . 20[ • We note however that the 3d wetting problem is special 
in that the critical behavior is non-universal. The critical exponents depend on the capillary 
parameter u [5j , which in turn involves the dependence on microscopic details of the system. 
In consequence, these cannot be viewed as completely irrelevant for the problem. 

One very convenient and widely applied cutoff choice, often referred to as Litim cutoff 



2lj is 

R A (q) = a(A 2 - q 2 )9(A 2 - q 2 ) , (17) 
which leads to the following flow equation 

d s u A (z) = - — g^T^x +{d~ l)u A {z) + ^-^zd z u A {z) , (18) 

where Ad = . This is not particularly surprising, that the non-linear part of the 

flow equation differs as compared to Eq. (fl4|) . More importantly however, upon expanding 
in d 2 u A the linear RG is again recovered, but the coefficient in front of the d 2 u A (z) term is 
now different by the factor B = -pj, yielding the conclusion that the capillary parameter 
oj ~ B^ 1 is a scheme-dependent quantity. 

This observation implies that the linearized RG framework involves a dependence of 
the crucial parameter of the theory on the approximation scheme. In the physically most 
relevant case d = 3 one finds B = 1 and consequently the two cutoff choices discussed above 
yield the same value of the capillary parameter. This is quite unlike the case d = 2, where 
the linear RG is however not an adequate approach. The present observation would give rise 
to sensitivity of critical behavior to the cutoff scheme for d g]2, 3[ when these are computed 
within the linearized theory. This dependence ceases as d — > 3~. For d > 3 the critical 
exponents carry no dependence on u are are determined by mean-field theory. 

The robustness of the capillary parameter in d = 3 against scheme variation can in fact 
be proven witout specifying the form of the cutoff function R A (q). To this aim we inspect 
Eq. f)12p and perform a linearization of the integrand 



o / s S d - 2 f , ±± d s r A (q) 
a s u A {z) = / dyy " 



4(2^)^7 l + r A (y) 



\ _ d 2 u A (z) 

(l + r A (y))y "\ 



3 d 

+(d-l)u A (z)^ — — zd z u A {z) . 



(19) 

The condition of scale invariance of the flow equation requires that r A (y) can be written as 
a function of exclusively one dimensionless variable (y), i.e. no additional dependence on A 
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occurs [lOj. We may write d s r^(q) = 2yd y r(y). The first integral on the RHS of Eq. ( TT9|) 
gives a physically irrelevant constant. The second one involves d^u\(z) multiplied by a 
cutoff-dependent constant, which determines the capillary parameter u. We immediately 



observe that for d = 3 



/ 



= 2r~' * =-2, (20) 
(1 + r A (y)) 2 ?/ A(^o) (1+r) 2 



where we used lim^o r(y) = oo and lim^oo r(y) = 0. This allows us to conclude that the 
scheme-dependence of the linear RG framework generically ceases as d — > 3~. 

The above fact and the observation that the terms neglected when performing the lin- 
earization of the RG flow equation involve only higher powers of the effective potential 
curvature, strongly support the reliability of the linear RG approach in d = 3. Cutoff de- 
pendence of the parameter u would give rise to a scheme dependence of the critical exponents 
when these are computed within the linear RG framework in d < 3. 

DISCUSSION 

We end by discussing the assumptions of locality and lack of renormalization of the stiff- 
ness coefficient a, which may seem controversial in view of more recent developments. In 
fact, the discrepancy between theoretical predictions based on the standard capillary- wave 
Hamiltonian Eq. ([T]) and numerical studies of wetting in the Ising model in d = 3 has 
led to significant refinement of the former. The history of the development is reviewed 
in Ref. [3|. An attempt to resolve the abovementioned discrepancy led to the conclusion 
that the capillary-wave Hamiltonian Eq. ([1]) should be amended by a weak dependence of 
the stiffness coefficient a on I jisj]. The RG studies of this model [lsj revealed that 
the transition is driven weakly first order for a range of parameters. However, this seem- 
ingly improved capillary-wave model strongly contradicted other firmly accepted aspects of 
wetting theory More recent developments suggest that a correct description of critical 
wetting phenomena requires a non-local interfacial model [if], \v^, which sheds light also 
on other puzzling features of wetting transitions including the wedge covariance property 



18 



19j and fulfillment of exact sum rules for the two-point correlation functions. The lin- 
ear RG analysis of the proposed non-local model supports the predictions of the classical 
capillary-wave model regarding the non-universal critical behavior, simultaneously revealing 
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the origin of the problem with the model retaining the dependence of a on I. Apart from 
subtleties, which however offer a solution to the long-standing problem why the Ising model 
simulations do not resolve the critical region, the local capillary-wave Hamiltonian Eq. (TjQ) 
turns out to serve as a good approximation to describe critical wetting behavior. 

The present work shows how the standard RG theory of wetting can be recovered relying 
on few clear assumptions. The observation, that the key parameter of the RG theory (u) is 
a scheme- dependent quantity makes the formal aspects of the approach somewhat weaker, 
but does not deprive it of predictive power in the vicinity of the physical dimensionality 
d = 3. Quite contrary, we showed that the ambiguities arise only below d = 3 and the linear 
RG generically becomes perfectly robust against scheme variation in the limit d — > 3~. 

I would like to thank Marek Napiorkowski for useful discussions and a number of com- 
ments on the manuscript. I also acknowledge a useful correspondence with Andrew O. 
Parry. 
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